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Abstract —In this paper, we study Two-Way Relaying (TWR) 
networks well-known for its throughput merits. In particular, we 
study the fundamental throughput delay trade-off in two-way 
relaying networks using opportunistic network coding (ONC). 
We characterize the optimal ONC policy that maximizes the 
aggregate network throughput subject to an average packet delay 
constraint. Towards this objective, first, we consider a pair of 
nodes communicating through a common relay and develop a two 
dimensional Markov chain model capturing the buffers’ length 
states at the two nodes. Second, we formulate an optimization 
problem for the delay-constrained optimal throughput. Exploit¬ 
ing the structure of the problem, it can be cast as a linear pro¬ 
gramming problem. Third, we compare the optimal policy to two 
baseline schemes and show its merits with respect to throughput, 
average delay and power consumption. The numerical results 
reveal interesting insights. First, the optimal policy significantly 
outperforms the first basefine with respect to throughput, delay 
and power consumption. Moreover, it outperforms the second 
baseline with respect to the average delay and power consumption 
for asymmetrical traffic arrival rates. 

I. Introduction 

Triggered by the seminal work of Alshwede et al. |l j, 
network coding has received considerable attention from the 
community, initially, to improve wired networks capacity j2j. 
Later, the merits of network coding prevailed in wireless 
networks due to the broadcast nature of the transmissions, 
opening up more opportunities for packet mixing |3)-(5). Re¬ 
cently, there has been increasing interest from the community 
in studying network coding in cooperative relay networks. 
Cooperation in wireless networks 0 constitutes one 
manifestation of spatial diversity which utilizes the broadcast 
nature of wireless transmissions to overcome the imperfections 
in wireless channels. TWR networks (8) have recently emerged 
as one of the basic forms of cooperative networks. 

Network coding schemes in TWR networks can be divided 
into two generic schemes referred to as 3-step and 2-step 
schemes 0. In 3-step schemes, each source node is allowed 
to transmit one packet to the relay in an exclusive time slot. 
In the third time slot, the relay broadcasts a bit-wise XOR- 
ed packet to both nodes to complete the two-way relaying 
process. On the other hand, 2-step schemes allow both source 
nodes to transmit simultaneously to the relay in one time 
slot. In the second time slot, the relay broadcasts the XOR- 
ed packet to both of them. In [9), the authors investigated 
and characterized conditions for maximizing the two-way rate 
for a number of 3-step and 2-step schemes. In addition. 



the work in 0 characterized the achievable rate regions 
for 3-step network coding in TWR networks. In 0, the 
concept of opportunistic network coding was first introduced 
to minimize the delay encountered by packets waiting at the 
relay node to be encoded. In p~2) , [13] the energy-delay trade¬ 
off was analyzed for conventional network coding. Most of the 
aforementioned works assume the presence of two buffers at 
the relay node in order to store packets from both sources 
to combine before relaying, using network coding. However, 
this setting leaves the multiple access channel (MAC) between 
the source nodes and the relay inefficiently utilized due to the 
stochastic nature of the packets arrivals at the source nodes. 
This, in turn, leaves no room for packet combining on the 
MAC channel resulting in throughput degradation. This work 
is motivated by the key observation that moving the buffers 
from the relay node to the source nodes would create more 
combining opportunities. This, in turn, results in “additional” 
throughput attained by other nodes in the network, as we show 
later in this paper. 

Our main contribution in this paper is three-fold. First, in 
contrast to 0, 0, we move the buffering storage from the 
relay node to the source nodes which allows us to enhance the 
throughput of the MAC channel, and hence, the entire TWR 
network throughput. Second, we characterize the optimal 
ONC policy, depending on the buffers’ states and the packet 
arrival rates at the source nodes, that maximizes the total 
network throughput (including the aforementioned additional 
throughput) subject to average packet delay constraints. Third, 
we investigate the average transmission power consumed at 
the source nodes for the proposed policy. Under the proposed 























setting for packet buffers at the source nodes, we show that 
the TWR network can sustain the same throughput achieved 
by conventional network coding, yet, at a lower average delay. 
Moreover, the MAC channel can be utilized more efficiently to 
communicate the traffic of other nodes in the network, giving 
rise to “additional” throughput. Towards solving the aforemen¬ 
tioned problem, we first construct a two-dimensional Markov 
chain capturing the buffers’ states at the two source nodes 
at hand. Afterwards, we formulate and solve an optimization 
problem that maximizes the total average throughput of the 
network under a constraint on the average end-to-end packet 
delay. Finally, given the established optimal TWR transmission 
policy, we characterize the minimum average transmission 
power consumed by the source nodes under this policy and 
compare it to the baseline policies. 

The rest of the paper is organized as follows. In Section 
II, the system model and underlying assumptions are intro¬ 
duced. In Section III, the opportunistic TWR transmission 
scheme is introduced. Afterwards, the constrained optimization 
problem is formulated and solved efficiently to characterize 
the maximum average throughput subject to average delay 
constraints,in Section IV. Also, the average transmission power 
for the established optimal policy is characterized in the same 
section. The numerical results are presented and discussed in 
Section V for a number of scenarios. Finally, our conclusions 
are summarized in Section VI. 

II. System Model 

We consider a TWR network consisting of two source 
nodes A and B communicating through a relay node R, as 
shown in Fig. [T] There is no direct link between node A and 
node B. Unlike prior work, each source node is assumed to 
have a limited size buffer; denoted Q a of maximum length 
N a at node A and Q\, of maximum length A), at node B. 
While the relay node, R, has no queuing capability. Thus, in 
essence, one of the byproducts of this work is to shed light 
and quantify the benefits of a simple design change, that is, 
moving the packet buffers from the relay node (as in ED- 
ED and many others) to the source nodes. As will be shown 
later, this simple design trick yields profound throughput gains 
attributed to a more efficient use of the MAC channel of TWR. 
In addition, we build upon this new system to characterize the 
optimal transmission policy at the source nodes to yield further 
performance gains beyond the state-of-the-art. 

We assume a time slotted system where each packet trans¬ 
mission fits exactly in one time slot. In addition, the two 
buffers’ backlogs are assumed to be available at the relay node 
which is assumed to take the transmission decision every time 
slot. The TWR communication process between nodes A and 
B is divided into two steps |5J; 

• Multiple access (MAC) step: nodes A and/or B transmit 
two/one packet(s) to the relay in one time slot. 

• Broadcast (BC) step: the relay R transmits the informa¬ 
tion back to the nodes in the subsequent slot. 

The packet arrival processes at source nodes A and B are 
independent and each follows an i.i.d. Bernoulli process with 


arrival rates A a and A& packets/slot, respectively, where 0 < 
A a , At, < 1. The channel between node A (or B) and relay 
R is assumed to suffer Rayleigh fading with scale parameter 
s a (or Sb). Let h ar denote the channel gain from node A to 
node R. Similarly, the channel gains from node B to node 
R, from node R to node A and from node R to node B are 
denoted by hb r , h ra and h r b, respectively. The channels are 
assumed to be reciprocal, i.e. h ar =h ra =h a and hb r =h r b=hb 
as shown in Fig. [T] The channel state information (CSI) is 
assumed to be fully known at the source nodes as well as 
the relay node. All nodes are assumed to be half-duplex, that 
is, a node can receive an arriving packet at the beginning of 
time slot, if any, and hold it in its queue, but cannot transmit 
a received packet within the same time slot. Let P a and I), 
denote the transmission power of nodes A and B, respectively. 
We assume that source nodes employ nested Lattice coding 
ED for simultaneous transmissions, such that in the event 
that two packets are received simultaneously at the relay node 
R , P a and Pb are chosen so that they can be readily decoded. 


III. Proposed Opportunistic network coding 
Two-Way Relaying Scheme 

In traditional network coding, e.g., ED’ ED’ the arriving 
packets are transmitted by source nodes A ( orB ) to node R, 
instantly, i.e. upon their arrival. In such setting, the packet 
buffers are assumed to reside only at the relay node, where a 
packet from one source always waits to be combined with 
a packet from the other source once it becomes available. 
However, this system setting has a fundamental limitation 
which causes the MAC channel to be busy at random time 
slots, governed by the packets arrival processes at nodes A 
and B, transmitting individual “uncombined” packets. This, 
in turn, presents a major source of inefficiency in the MAC 
channel which we successfully remedy in this paper by using 
2-step rather than the 3-step network coding approach. 

The proposed scheme makes more efficient use of the 
MAC channel, since the buffers are assumed to be at the 
source nodes (instead of the relay node). This, in turn, makes 
it possible to control the time slots in which the MAC 
channel is occupied. This work quantifies the performance 
gains attributed to the optimal TWR transmission scheme 
under the new setting, in the sense of maximizing the total 
network throughput under an average delay constraint. We also 
characterize the average transmission power consumed by the 
optimal policy compared to two baseline schemes, including 
traditional network coding 1121. 


A. Queuing-theoretic Model 

In order to analyze the performance of the proposed 
scheme, the buffers at nodes A and B are modeled as a 
two-dimensional Markov chain with steady state probabilities 
iTij where i and j represent the number of packets queued 
at nodes A and B, respectively. The states of the Markov 

chain are denoted by S(i,j), where i = 0,1,2,. N a and 

j = 0,1, 2,. Nb- The proposed randomized policy, assumed 

to be carried out by the relay node, takes state-dependent 





transmission decisions. Define of-, k = 1,2,3,4, to denote the 
transmission probabilities under the four cases of interest: 

1) Only node A transmits one packet: with probability g]-. 

2) Only node B transmits one packet: with probability gf.-. 

3) Nodes A and B transmit two packets simultaneously 
using Lattice coding with probability gfj. 

4) Neither node A or B transmit with probability gfj. 
with the constraint that these probabilities should sum to one, 
that is, for each state S(i,j), 
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E4 = 1 

it*—1 

(1) 

Next, we characterize the probabilities of packets arrivals at 
the two source nodes, A and B, in an arbitrary time slot. This 
gives rise to one of four cases, namely one packet arrival at 

each source, only one packet arrival at node A, 

one packet 

arrival at node B and, finally, no packet arrivals. 

Hence, we 

define /,;, i = 1,2,3,4, to denote these four events. 

in order: 

fl = X a X b 

(2) 

,f'2 = A a (1 - X b ) 

(3) 

,/*3 = (1 — Aa) A b 

(4) 

h = (1 - Aa) (1 - A 6 ) 

(5) 

Next, we characterize the state transition probabilities of the 
Markov chain at hand using the characterized arrival rates at 
Q a and Q b in (2)-(5), along with the transmission probabilities 

gf. as follows. 


First, if both buffers are empty (i.e. the origin state), it is 
straightforward to notice that we have four possible transitions. 

P(5(0,0)|5(0,0)) = / 4 

(6) 

P (5 (0,1)15(0,0)) = A 

(7) 

P (5 (1,0)|5(0,0)) =/ 2 

(8) 

P (5 (1,1)|5(0,0)) = fi 

(9) 

Second, if only Q a is non-empty which corresponds to states 
on the horizontal axis, S ( i, 0), * = 1,2, 3, . N a , 

P(S(i,0)\S(i,0)) = f 2 gl o + Ugf o 

GO) 

P (5(t + 1,0) \S(z,0)) = f 2 gf o 

(11) 

P(S(i-l,0)\S(i,0)) = f 4 gf 0 

(12) 

P(S(i,l)\S(i,0)) = figl 0 + f 3 gf 0 

(13) 

P {S (i + 1,1) 5 (i, 0)) = figf 0 

(14) 

P (S {i ~ 1,1) 5 (i, 0)) = f 3 g} 0 

(15) 

As for the vertical axis states, S(0,j), j = 1,2, 

3,. N b , 


the state transition probabilities can be derived along the 
same lines of (10)-(15). Finally, we consider the general case 
where both node A and B have non-empty buffers, i.e. the 

interior of the state space, S ( i,j ) for i = 1,2, 3,. N a ,j = 

1,2,3,. N b , 

P{S ( i,j) | S (■ i,j)) = f 4 gtj + figfj + hg% + hgfj (16) 
P(S(i~l,j-l)\S(i,j)) = f 4 gfj (17) 

P(S{i + l,j + l)\S(i,j)) = figfj (18) 
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Given the introduced Markov chain model, we formulate and 
solve the target optimization problem in the next section, 
aiming at maximizing the total network throughput under an 
average packet delay constraint. 

IV. Optimal Throughput Policy under Delay 
Constraints 
A. Problem Formulation 

In order to demonstrate the throughput gains on the MAC 
channel in our setting, which contributes to enhancing the total 
TWR throughput, we assume the presence of another source- 
destination pair denoted by CD, respectively. The pair CD is 
assumed to utilize the MAC channel idle slots, i.e. not used 
by pair AB. These nodes may be readily thought of as low 
priority opportunistic users who are continuously on the look 
for spectrum holes (unused slots in our context), e.g. cognitive 
radio users with perfect sensing capability. 

Now let fi \ denote the throughput of nodes A and B while 
fi 2 denotes that of nodes C and I), 

Pi = ^ * ((sh + 9ij) + 2 * 9ij) ( 25 ) 

Pi='52 7 Tij * gfj (26) 

i,j 

Ptot = Pi + P-2 (27) 

We assume buffers Q a and Q b to be non-lossy thus (28) will 
be always satisfied. 

Pi = A a + Ah (28) 

Thus by moving the buffers from the relay node to the source 
nodes A and B, we can control the time slots in which the 
MAC channel is empty along with maintaining (28). 

It is evident now from the above discussion that maximizing 
the total network throughput is equivalent to maximizing the 
number of empty slots on the MAC channel, represented by 
7 tij gfj- This, in turn, gives rise to an instance of the funda¬ 
mental throughput-delay trade-off which recurs in different 
problem contexts in wireless communications and network¬ 
ing. Motivated by this trade-off, we formulate a constrained 
optimization problem for maximizing the number of empty 
slots (equivalent to total throughput as shown in the previous 
discussion) subject to a constraint on the average packet delay, 
as follows: 

PI max Tijgfj (29) 

ZJ 1,3 

S.t. - —V' 7 Tij (■ i+j)< D max (30) 

A a i Afr 

i,3 







TrP = 7T, ^ Kij = 1 

(31) 

l i3 

9oj = 9io = 9oj = 9io = 0 

(32) 

9% a j = 9iN b = 1 

(33) 


where D max is the average packet delay constraint and P is 
the state transition probability matrix of the system Markov 
chain described in Section III. It should be noted here that 
the objective function in PI is the average number of empty 
slots, characterized by the no transmission probability gfj 
introduced earlier. The first constraint is the average packet 
delay constraint obtained using Little’s law. The constraints 
in ( |Ti~[ ) are the balance equation and probability normalization 
condition for the Markov chain model, respectively. Then, in 
order to govern the physics of the system, the probabilities in 
([32} are set to zero indicating no transmissions are possible 
from empty buffers. Also, to avoid any packet loss at the 
buffers, a source node always transmits a packet from a full 
buffer with probability one as in ( [33} . 

B. Problem Complexity and Solution Approach 


In order to solve the above problem we need to obtain the 
steady state probabilities of the Markov chain governing the 
system dynamics, 7r, ? . Motivated by the sheer complexity of 
getting a closed form expression for the system Markov chain 
to be plugged in PI, we resort to a change of variables similar 
to G3 in an attempt to linearize the problem and, hence, 
significantly simplify the solution. 

To this end, we introduce the intermediate variables x k 3 

where for i = 1,2,3. N a , j = 

1,2,3,. Ni> and k = 1,2,3,4. Thus, we can solve the 

transformed problem for and then map it back to the 
transmission probabilities, g^, using (34 1 . 


9ij = 


(34) 


V r k 

2-^k—l x ij 

Accordingly, problem PI can be transformed to the equivalent 
problem P2 as follows. 


P2 

max > x% 

x k - ^ 3 

(35) 

s.t. 

x a + x b ^^ Xi ^ l+J )- Dmax 

(36) 


i,j k—1 

4 

Q x = 0 ’ = 1 

(37) 


i,j k—1 

X M = x i0 = x 0j = x 0j = 0 

(38) 


X Naj = X N a j = X Naj = ^ 

(39) 


x iN b = x lN b = x iN b = 0 

(40) 


where x is the vector of the new optimization variables 

[ x oo, x oi, x oi, x oi> x on . ’ x ij’ x ij’ x ijy x iji .] and 

Q is the transition equations matrix in terms of the new 
variables a :')). In addition, we have the boundary conditions in 
(|38}-([40} corresponding to the boundary conditions (32 1 ,(331 
of problem PI. 


321,(33 


C. Average power consumption for the optimal transmission 
policy 


Given the optimal transmission policy characterized in Sec¬ 
tion IV-B to maximize the total network throughput, our prime 
goal in this subsection is to determine the average transmission 
power for each node depending on the available CSI. In 
order to sustain a constant transmission rate at the source 
nodes, the consumed power varies according to the varying 
channel gain, i.e. a node should increase its transmission 
power during poor channel conditions. Thus in order to avoid 
wasting high transmission power, we propose that a packet can 
only be transmitted from a node if the channel gain for the 
corresponding link is greater than or equal to a threshold hth- 
For the assumed Rayleigh distributed channel gain realizations 
with unit scale parameter s = 1, the probability that the 
channel gain h is greater than or equal to h t h in a given time 
slot is given by: 


P(h> hth ) = exp - 


h 2 th 


(41) 


Note that according to the theory of Lattice coding fl4) , it has 
been established that the source nodes can effectively achieve 
a rate of | log (0.5 + h* P) for a channel gain of value h 
and transmission power of value P. Assuming normalized 
noise power at the source nodes A and B, let a = (2 2r — l) 
and /3 = (2 2r — 0.5) denote the SNR in case of a single 
transmission or simultaneous transmissions, respectively. And 
r represents the rate of transmission per time slot, where r is 
1 when a packet is transmitted and zero when no packets are 
transmitted. 

Recall that the optimal transmission probabilities g k * result¬ 
ing from P2 are state dependent. Thus, in order to minimize 
the average transmission power, we map these probabilities to 
a channel threshold for nodes A and B as follows: 



Note that (|42}-(|45), can be derived for a Rayleigh distributed 
channel, i.e. gj*=P [h a . i3 > h tha i .).P (h btij < h thbi .)\/ i,j. 
Averaging over all possible realizations of the Rayleigh dis¬ 
tributed channel, the average power consumed at the source 
nodes in every time slot follows one of the three possibilities 
i) If h a ,ij (hb^ij') hth a (b),ij &: h bl ij {ha,ij')^h / th b ( a '),ij then. 


Pa,ij{Pb,ij) — a \/ 2 61 ^ C 


hth r 


( W >«1 


V2 


1- exp - 


h 2 


Pb,ij (Pa,ij) = 0 
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Fig. 2: Total average throughput of nodes A, B, C and D, A a =A&. 



Fig. 3: Total average delay at nodes A and B, A a =A b - 


ii) If ha,ij ^ ^th a ,ij hb,ij 'Hh'thb,ij then, 


R. The first baseline scheme denoted by Random MA is 


p a^=P ■ / |erfc(^j exp - 2 


Pb.ij = P\ 77 erfc 



h 2 

rL th b ,ij 


exp 


h 2 

rt th a ,ij 


iii) If da,ij < ^'h j th a ,ij db,ij < ^-h,tht ) ,ij then, 


Ph.ii 0 


the same as that analyzed in 1121, where node R is assumed 


(48) 

(49) 

(50) 


Using the above transmission policy, would significantly 
decrease the transmission power consumed by the source 
nodes as it aims to prevent wasting unnecessarily increased 
transmission power to overcome bad channel conditions. On 
the other hand, although the simultaneous transmissions would 


slightly increase the transmission power 1141, the proposed 
policy results in lower average transmission power compared 
to conventional network coding as shown in section V. 

In the next section, we present our numerical results show¬ 
ing the merits of the optimal TWR transmission policy with 
respect to maximizing the total network throughput subject to 
average delay constraints in addition to the average transmis¬ 
sion power consumption by the source nodes. 

V. Numerical Results 

The numerical results presented in this section summarizes 
our major findings and lessons learned from this work. It is 
straightforward to establish that the formulated optimization 
problem P2 is a linear program that can be solved efficiently. 
In order to show the merits of the optimal policy, we will 
compare its performance to two baseline schemes, namely 
conventional network coding 03. fl3| , in terms of the total 
average throughput, average packet delay and the average 
transmission power consumption at the source nodes. 

As for the conventional network coding baselines, we 
assume having two systems, each consists of two source 
nodes A and B exchanging information through performing 
conventional network coding at an intermediate relay node 


to have two finite length buffers. The source nodes transmit 
the arriving packets, whenever available, to node R, then 
every packet arriving at node R waits to be combined with 
another packet from the opposite traffic direction. While in 
the second baseline denoted by Combined MA, we assume 
having the same model but with the buffers now present in 
the source nodes instead of the relay thus the randomness in 
the MAC phase is no longer present. Note that in all of the 
three schemes, the broadcast channel operation is the same, 
such that the relay node broadcasts the received packet(s) after 
performing the Lattice coding operation if required. 

We assume similar parameters in all the compared schemes 
as follows; the buffers at both source nodes are of equal 
length where N a = N h = 15 packets, packets arriving at 
nodes A and B follow an i.i.d Bernoulli process with equal 
average arrival rates A a = Af, packets/time slot. The total 
delay constraint applied to the proposed transmission scheme 
Dmax = 3 time slots. While in conventional network coding, 
the arriving packets from both directions are always queued 
awaiting a combining opportunity whenever available which 
induces higher delay. Since all three transmission schemes 
studied are inherently non-lossy, the average sum throughput 
of the two source nodes A and B is equal to the sum average 
arrival rates. 

In order to show the benefit attributed to increasing the 
average number of empty slots in the MAC channel by moving 
the buffers from the relay node to the source nodes, we will 
consider the presence of the additional nodes C and D, with 
perfect sensing capabilities as described previously in section 
IV. These nodes try to access the same channel only during 
these empty slots, to exchange some information packets of 
their own. For fair comparison, we assume that nodes C and 
D are also present in the other two baseline systems Random 
MA and Combined MA and operate only during empty slots, 
in a manner exactly similar to our proposed policy. 

Fig. [2] shows the total average throughput fi tot of the overall 
system (four nodes: original nodes A and B along with the 
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combine all the arriving packets before transmission. 



Fig. 4: Average power of nodes A and B, Ab=0.5. 


opportunistic users C and D ) making use of the spared empty 
slots in the given system. We assume that nodes C and D 
always have queued packets to transmit. It can be easily no¬ 
ticed that the total network throughput significantly increases 
when the buffers are present at the source nodes where the 
MAC channel can be efficiently used. The total throughput 
achieved from the proposed optimal scheme approaches this 
of the second baseline. Combined MA. In this case, nodes A 
and B can transmit simultaneously, large portion of the time, 
which also increases the number of empty slots on the average 
for node C and D. 

In the baseline systems using Random MA and Combined 
MA, the queued packets incur the same delay as they both 
use conventional network coding whether at the node R or at 
nodes A and B. And since the proposed transmission scheme 
is designed under a total delay constraint, Fig.[3]shows that the 
previous throughput gain can be acquired under the required 
delay constraint. Thus our transmission scheme approaches the 
total throughput gain of system using Combined MA, giving 
the maximum throughput gain, however with significantly 
lower delay. 

The second part of the results is about the average power 
consumed at nodes A and B for asymmetrical arrival rates, 
\ a = 0.5 and A b € [0.1,1] packets/time slot. Given the 
transmission probabilities of the proposed scheme, the mini¬ 
mum average power consumption at nodes A and B is shown 
in Fig. [4] based on the mapping mentioned previously in 
section IV. Since the second baseline. Combined MA, uses 
conventional network coding and the source nodes are aware 
of simultaneous transmission, the average power consumed is 
maximum with different arrival rates. On the other hand, the 
source nodes in Random MA transmits their packets to the 
relay node upon arrival since the relay is the node responsible 
for combination. Thus nodes A and B do not have to increase 
their transmission power in this case consuming lower power 
consumption. Finally, the proposed ONC scheme achieves the 
least power consumption levels for asymmetrical arrival rates 
approaching that of Combined MA only at equal arrival rates, 
A a = Af, = 0.5 packets/time slot where the system tends to 


VI. Conclusion 

In this paper, we studied the throughput delay trade-off 
in TWR networks using ONC. Specifically, we considered a 
pair of nodes, each having a finite length buffer, communi¬ 
cating through a common relay node. We characterized the 
optimal ONC transmission policy that maximizes the aggre¬ 
gate network throughput subject to an average packet delay 
constraint. First, we formulated an optimization problem for 
the delay-constrained optimal throughput. Then, we compared 
the performance of the proposed transmission policy to two 
other baseline schemes. The numerical results showed that the 
optimal policy outperforms the first baseline with respect to 
the average throughput, delay and power consumption. Also, 
it outperforms the second base line with respect to the average 
delay and power consumption for asymmetrical traffic rates. 
It is of interest as a future work, to investigate the general 
optimal transmission policy where the source nodes as well 
as the relay node both have queuing capabilities. 
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